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We present a new open-source software package, FIT-MART, that allows non-experts to quickly 
get started simulating quantum magnetism. FIT-MART can be downloaded as a platform- 
idependent executable Java (JAR) file. It allows the user to define (Heisenberg) Hamiltonians 
by electronically drawing pictures that represent quantum spins and operators. Sliders are auto- 
matically generated to control the values of the parameters in the model, and when the values 
change, several plots are updated in real time to display both the resulting energy spectra and the 
equilibruim magnetic properties. Several experimental data sets for real magnetic molecules are 
included in FIT-MART to allow easy comparison between simulated and experimental data, and 
FIT-MART users can also import their own data for analysis and compare the goodness of fit for 
different models. 
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I. MOTIVATION 

For thousands of years, people have known that certain materials are magnetic, but the origin of this phenomenon 
remained a mystery until the last century. This is because any description of material magnetism relies on concepts 
from the (modern) fields of quantum and statistical physics. Today, this phenomenon is well-understood by condensed 
matter physicists; but it still remains a mystery to most people, including most undergraduate physics students. This 
is the problem that we seek to address. To make quantum magnetism accessible to a broader audience, we have created 
a user-friendly, open-source software package for modeling simple models of quantum magnetism. This software can 
be downloaded from the Open Source Physics collection pQ as a platform-independent executable Java (.JAR) file. To 
get started, a user simply downloads this file and double clicks. 

The software introduced here will be referred to with the acronym FIT-MART which stands for "Fully-Integrated 
Tool for Magnetic Analysis in Research and Teaching". It is a "fully integrated tool" in the sense that the entire 
modeling process is self-contained within one program. There is no need to first compile the code; results are 
immediately provided in the form of several plots; and these data can be directly compared to experimental data — all 
within a single window on a computer. (Results can also easily be exported as images or text files for later use.) 
This software is also sufficiently sophisticated that it can provide non-trivial results for certain types of real magnetic 
molecules, and it is currently used by several researchers who study molecular magnetism. To complete the acronym, 
FIT-MART has been designed to be a viable instructional tool for undergraduate courses (e.g., modern physics, 
statistical mechanics, and condensed matter physics), even for students who do not have previous computational 
experience. This is demonstrated for a few specific examples in Sec. |III| 

It is important to note that many advanced simulations are available for researchers who study condensed matter 
systems. In the field of magnetism, the ALPS project [2^4] provides very powerful codes and libraries for simulating a 
wide variety of quantum lattice models (Heisenberg, Hubbard, etc.) using several different numerical methods (exact 
diagonalization, quantum Monte Carlo, Density Matrix Renormalization Group, etc.). This versatility of models and 
methods gives ALPS a major advantage compared to FIT-MART, but this versatility comes at a price: FIT-MART 
is much simpler to use. The strengths and limitations of FIT-MART are explored below in Sec. [TTJ 

II. IMPLEMENTATION 

In order to achieve simplicity in FIT-MART — yet still provide non-trivial results — all simulations use models that 
are described by the Heisenberg Hamiltonian, 

(hj) i 

where tildes denote quantum operators. Here Sj and Sj represent the dimensionless spin operators [5] corresponding 

to two sites, i and j, and Jij describes the strength of the interaction (bond) between these two sites (J > for an 
antiferromagnetic interaction). The second term gives rise to the Zeeman effect: g is the spectroscopic splitting factor, 
Hb is the Bohr-magneton, and H is the external magnetic field. Defining H to point along the z-axis, the second 
term becomes —gjiBHS Zl where S z = ^ and Siz represents the z component of the individual spin operator for 

site i. This Hamiltonian commutes with both of the total spin operators, 5 2 (where S = Sj) and S z , which have 
the familiar eigenvalues 5(5 + 1) and Ms, respectively. To produce results, we generate the matrix representation of 
Eq. (fil using the basis wherein the individual spin operators, s^, are diagonal, and then exact diagonalization is used 
to calculate both the energy eigenvalues and eigenvectors. Since S z commutues with we are able to separately 
treat each subspace of % that has a different Ms eigenvalue. This reduces the size of the matrices to be diagonalized, 
but FIT-MART is still limited to relatively small numbers of spins. 

To demonstrate the practical limations of FIT-MART for current personal computers, we have measured computa- 
tion times, t, for several different system sizes. [6] These times are shown in Table [l[ where N represents the number 
of spins, each of spin s. The dimension of the total Hilbert space is then D = (25 + 1)^, and D 1 is the dimension of 
the largest subspace — i.e., the largest matrix to be diagonalized is D' x D' . When t <C 1 sec, there is no perceptible 
time delay between the user dragging a slider (e.g., to update Jij) and the plots updating. Hence, FIT-MART runs 
smoothly for t <C 1 sec, and there is a noticeable delay for larger system sizes. For D' > 1000, the calculation of all of 
the energy eigenvalues requires more than 2 GB of RAM, so these computations become impossible on computers that 
do not have sufficient memory available. To see exactly how FIT-MART is implemented, we encourage the interested 
reader to download the source code from Ref. pQ. FIT-MART was created using Easy Java Simulations which is a 
Java code generator that is popular in computational physics education. [7H9] 
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TABLE I: Computation times for different system sizes. [6] 
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III. MODELING 



We now demonstrate how FIT-MART is actually used. Four screenshots from FIT-MART are displayed in Fig. [T] 
Fig. [TJa) shows a user-drawn picture that corresponds to the Hamiltonian 



% Jisi • s 2 + J 2 (si • s 3 + 5 2 • s 3 ) - gjJLBH ■ (s lz 



S2z 



(2) 



where the top (red) circle represents S3 and the bottom (blue) circles represent si and 82. A user can draw a 
Hamiltonian such as this by first clicking the buttons shown in Fig. [ljb) — e.g., "Add Spin" and "Add Bond" — and 
then using a mouse to properly position the spins and bonds within Fig. [IJa). After defining the structure of the 
Hamiltonian, clicking "Calculate and Plot" will produce the controls that are shown in Fig. [ljc), which define the 
values of the parameters that appear in Eq. [2j and it will also produce the controls shown in Fig. [TJd), which define 
the temperatures and fields that are used to calculate equilibrium magnetic properties. 




(a) 



Bonds: | Kelvin S ] 



Jl = 10.00 Kelvin 



Min = -25.00 Kelvin 



Max = 25.00 Kelvin 



JL = L0.00 Kelvin (AFM)- 



-J2 = -5.00 Kelvin (FM)- 



2.00 



g = 2.00 



O 



^Type A spins: s = 2.0- 



-Type E spins: s = 2.5- 



=0= 



0.5 1.5 2.5 3.5 4.5 



(c) 



Remove 5 pin | | Add Spin | of type: [ A ± | 



Remove Bond | | Add Bond | of type: [ Jl 



Save Molecule 



Load Molecule 



Clear Molecule 



HELP! 



Calculate & Plot 



Load All 



Mag Parameters 



Temp = 0.50 K 



Min: 0.10 K 



Max: |4.00 | K 
Temp = 0.5 Kelvin 



o 



Max Field: 20.00 T 



Min Field: 0.00 T 



Change in Field: 0.05 T 



X Parameters 



Field = 0.10 T 



Min: 0.00 T 



Max: 7.00 T 



Field = 0.1 Tesla 



o 



Max Temp: 1 300.0 | K 



Min Temp: |0.10 | K 



Change in Temp: 0.10 ^K 



(b) (d) 
FIG. 1: Screenshots from FIT-MART showing how a user defines a model to be simulated. 
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Field (Tesla) Field (Tesla) Field (Tesla) 

(a) (b) (c) 

FIG. 2: FIT-MART output for a dimer of spins s = 1/2 with J/k B = 15 K, and g = 2. For (b) and (c), curves are 
shown for three temperatures: 0.5 K (blue), 2 K (red), and 4 K (green). 



Fig. [2] shows FIT-MART output for an extremely simple example: two spins, each with s = 1/2, coupled by 
an antiferromagnetic bond (J > 0). This example is pedagogically useful for understanding several topics from 
introductory quantum and statistical physics. In Fig. [2|a), we can see an S = singlet (the zero-field ground state) 
and an S = 1 triplet, where the zero-field energy gap [10 between the singlet and the triplet is proportional to 
J. The zero-field degeneracy of the S = 1 triplet is lifted by the external field with a slope (of E vs. H) that is 
proportional to g, and when H « 11 Tesla, the ground state changes between \S = 0, Ms = 0) and \S = 1, Ms = 1). 
The effect of this ground state "level crossing" is clearly evident in the plots of magnetization, M, and differential 
susceptibility, dM/dH, versus field, H. There is an abrupt step in M(H) at low temperatures, T, and a corresponding 
narrow peak in dM/dH; with increasing T, this step (and corresponding peak) are broadened. These features are 
entirely determined by the Maxwell-Boltzmann distribution function: At very low temperatures, only the ground 
state is thermally populated; and at elevated temperatures, excited states become increasingly populated. Hence, this 
example provides a simple platform for student exploration in statistical physics. 

In addition to the data that are plotted versus H in Fig. |2j FIT-MART also calculates and plots three quantities 
versus temperature for fixed (small) values of H: x = M/H, Tx, and 1/x (not shown here). For this particular model 
(two spins s = 1/2 with antiferromagnetic coupling), the plot of x(T) has a peak at T « | J/kB- A peak in x(T) is a 
typical signature of antiferromagnetism, so it is desirable for students to understand the origin of this peak; and this 
is easily facilitated with FIT-MART. From Fig.[2|b), the slope of M(H) — which is x — has a value of x ~ for H < 5 
Tesla and T < 2 Kelvin, and x increases with increasing T. If you continue to increase T [using the slider labeled 
"Temp" in Fig. [TJd)], you will see that the slope of M(H) continues to increase until T ^ 9.3 K and then decreases 
as T is increased further. These results define the plot of x(T), which has a peak at T w 9.3 K. The plot of Tx vs. T 
increases monotonically with increasing T, approaching the Curie constant, C, for two independent spin 1/2 particles 
as T 00. In fact, T X > 0.99C for T = 300 K. 

Fig. [3] shows FIT-MART output for a less simple, experimentally relevant example. Fig. [3Va) shows published [TT] 
Tx data (symbols) for a {Mn3} molecule whose Hamiltonian is defined by Fig. [TJa) and Eq. Using the notation 
of Eq. si = S2 = 2 (there are two Mn 3+ ions), and 53 = 5/2 (there is one Mn 2+ ion). The best fit for these 
data [the solid curve in Fig. [3^a)] is achieved using the parameters Ji/ks = 9.58 K and J2/&B = 17.04 K, and a 
quantitative measure of the goodness of fit (the root-mean-square deviation between experiment and theory) for this 
set of parameters is displayed within Fig. [3^a). FIT-MART also allows one to include corrections for the bulk materal 
including diamagnetism and impuritities^l2 and diamagnetism has been included as a fitting parameter to obtain 
the data in Fig. [3^a). This set of parameters produces the zero-field energy spectrum shown in Fig. [3^b), which has 
an S = 3/2 ground state. Fig.[3^c) contains predictions for high-field, low-temperature magnetization measurements. 
If this model is correct, five ground state level crossings should be observable between H = 30 and 80 Tesla, resulting 
in the eventual saturation of the magnetization (to an S = 13/2 ground state) when H « 78 Tesla. We note that 
these are new predictions that could be tested with pulsed-field measurements. 
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1 2 3 4 5 6 7 

Temperature (Kelvin) Total Spin Field (Tesla) 

(a) (b) (c) 

FIG. 3: FIT-MART output for the Hamiltonian defined by Fig.Qa) and Eq. Q for J Y jk B = 9.58 K, J 2 /k B = 17.04 
K, and g = 2.1. (a) Measured [TTj (symbols) and calculated (curve) T%(T) data for a fixed field, H = 7 Tesla. A 
diamagnetic correction was included 4 as an additional fitting parameter, which yielded the value Xd = 1-5 x 10 -3 
cm 3 /mol. (b) The zero-field energy spectrum versus the total spin quantum number, 5, for the model that provides 
the best fit for T\- (c) dM/dH versus H for three temperatures: 0.5 K (red), 2.0 K (green), and 5.0 K (blue). 



IV. SUMMARY AND FUTURE DEVELOPMENT 



We have introduced an easy-to-use software package ( "FIT-MART" ) that allows users to perform advanced quantum 
mechanics calculations and analyze magnetic data for certain classes of real materials. This software is sufficiently 
user-friendly that it can be used in undergraduate instruction without previous computational experience. It is also 
sufficiently powerful that it has been used by several researchers studying magnetic molecules. An important feature 
of FIT-MART is the inclusion of experimental data, and the real-time analysis of the goodness of fit. [See Fig. [3ja).] 
The process of varying model parameters and comparing simulated and experimental data is an important part of 
scientific discovery, and FIT-MART provides a means for undergraduate physics students to be able to experience 
this modeling process first hand. [13] 

In subsequent versions of FIT-MART, the calculation of energy eigenvalues will be parallelized, which will sig- 
nificantly reduce the computation times listed in Table [I] when run on a multi-core computer. Automated fitting 
routines will also be incorporated in order to automatically search multidimensional parameter spaces. Finally, a 
series of exercises will be included within the FIT-MART documentation in order to support instructors wishing to 
use FIT-MART within their courses. 
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